Spurious detection of phase synchronization in coupled nonlinear oscillators 
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Coupled nonlinear systems under certain conditions exhibit phase synchronization, which may change for 
different frequency bands or with presence of additive system noise. In both cases, Fourier filtering is tradition- 
ally used to preprocess data. We investigate to what extent the phase synchronization of two coupled Rossler 
oscillators depends on (1) the broadness of their power spectrum, (2) the width of the band-pass filter, and (3) the 
level of added noise. We find that for identical coupling strengths, oscillators with broader power spectra exhibit 
weaker synchronization. Further, we find that within a broad band width range, band-pass filtering reduces the 
effect of noise but can lead to a spurious increase in the degree of synchronization with narrowing band width, 
even when the coupling between the two oscillators remains the same. 
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In recent years both theoretical and experimental studies of 
coupled nonlinear oscillators has demonstrated that such os- 
cillators can exhibit phase synchronization [1-5]. Analysis 
of experimental data has also indicated the presence of phase 
synchronization in a range of coupled physical, biological and 
physiological systems [6-17]. In many of these studies, an 
important practical question is how multi-variate time series 
characterized by relatively broad power spectrum are phase 
synchronized in a specific frequency range [18-24]. The pres- 
ence of internal or external noise may also be an obstacle 
when quantifying phase synchronization from experimental 
data [18, 19, 25-27]. In both cases a band-pass filter is tra- 
ditionally applied either to reduce the noise effect or to extract 
the frequency range of interest. Thus, it is important to know 
to what extent the width of the band-pass filter influences the 
results of the phase synchronization analysis, as well as what 
is the range of the index values obtained from the analysis that 
indicate a statistically significant phase synchronization. 

To address these questions, we consider a system of two 
coupled Rossler oscillators (1,2) defined as 



a:^i,2 = -t^i,22/i,2 - ^1,2 4 
2/1,2 = ti;i,2a;i,2 + ayi,2, 
z\,i = ./ + zi,2(a;i,2 -h) 
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with parameter values a = 0.165, / = 0.2, and h ~ 10. 
For the mismatch of natural frequencies, we choose cji 2 = 
Wo ± Atj, with Wo = 0.6 and Aw = 0.005 [Fig. 1(a)]. The 
time step in our simulation is Ai = 27r/10'^, and the signal 
length n = m<[t/At] with t = 10^, where int[x\ denotes the 
integer part of x. 

We first investigate the characteristics of the system defined 
in Eq. (1) by comparing them with the characteristics of a 
second set of two coupled Rossler oscillators (3,4) studied in 
[3]. The system (3,4) is also described by Eq.(l), and has the 
same values for parameters a, /, and h as system (1,2). The 
only differences are the natural frequency wo = 1 and the fre- 
quency mismatch Aw = 0.015 [Fig. 1(b)]. We observe a sig- 
nificantly broader power spectrum for system (1,2) with wo = 



0.6 and frequency mismatch Aw = 0.005 [Fig. 1(c)]. Further, 
we observe that the instantaneous phase differences A-0i 1 = 
[<t>xi {t) — <l>x2{t)] mod (27r) for system (1,2) exhibits larger 
fluctuations [Fig. 1(d)], described by a broader distribution 
[Fig. 1(e)], compared to system (3,4), suggesting a weaker 1:1 
phase synchronization for system (1,2). To quantify the de- 
gree of phase synchronization in the two Rossler systems we 
use the synchronization index p = (S'max — <S')/>5'inax [18], 
where S = — J2k=i ^^ InP/s is the Shannon entropy [28] of 
the distribution P(AV'i,i) of AVji,!, and S'max = In TV, where 
N = Mii[cxp(0.626 + 0.41n(n— 1.0))] is the optimized num- 
ber of bins over which the distribution is obtained [29]. For 
system (3,4) with a narrow power spectrum we obtain a sig- 
nificantly larger value of p compared to the system (1,2) char- 
acterized by a broader power spectrum [Fig. 1(f)]. Varying the 
values of the coupling strength C, we find that the phase syn- 
chronization index p is consistently higher for system (3,4) 
characterized by the narrower power spectrum. Thus, for the 
same coupling strength C and for identical other parameters, 
system (1,2) with wo = 0.6, which has a broader power spec- 
trum, exhibits weaker synchronization compared to system 
(3,4) with Wo = 1, which has a narrow power spectrum. These 
findings are complementary to a recent study indicating a dif- 
ferent degree of phase synchronization for the spectral com- 
ponents of coupled chaotic oscillators [30]. 

Recent work has shown that coupled Rossler oscillators 
may exhibit different degrees of synchronization for differ- 
ent ranges of time scales obtained via wavelet transform [31]. 
Here, we ask to what extent the width of a band-pass filter 
affects the degree of phase synchronization between two cou- 
pled Rossler oscillators. While the output observables xiand 
X2 of system (1,2) are clearly not in phase [Fig. 2(a)], after 
Fourierband-pass filtering in the range of A/ = 0.01 centered 
at the peak of the power spectrum 27r/ « 0.54 [Fig. 1(c)], 
the observables xi and X2 appear 1:1 synchronized with well 
aligned peaks [Fig. 2(b)]. The effect of the band-pass filter 
can be clearly seen in the behavior of the instantaneous phase 
difference Ai/)i:i [Fig. 2(c)] and in the shape of the probabil- 
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FIG. 1 : Differences in tiie synchronization of two Rossler systems with identical coupling strengths and different power spectra. Phase plot 
trajectories of the variables x vs. their Hilbert transform x for: (a) system (1,2), with xi corresponding to uji = wo + Acu, where cua = 0.6 
and Alo — 0.005; (b) system (3,4), with xs corresponding to ua = luq + Aui, where now a;o = 1 and Aco — 0.015. For both Rossler 
systems C = 0.03. (c) Power spectra of the time sequence xi (dashed line) and xa (solid line). A broader spectrum is observed for system 
(1,2) compared to system (3,4). (d) Instantaneous phase difference A?/)i,i = {<pxi{t) — </'i,(i)) mod (27r) for system (1,2) (dashed line), 
and A'i/'i,! = {4>x3{t) ~ 4>xi(t)) mod (27r) for system (3,4) (solid line), and (e) their corresponding distributions P(A'i/)i,i). System (1,2) 
exhibits larger fluctuations in A'i/'i. i and is characterized by a broader distribution P{Ail}\.i). (f) Synchronization index p as a function of 
the coupling strength C. For identical values of C, system (3,4) (solid line) which is characterized by a narrower power spectrum exhibits 
stronger synchronization (larger index p) compared to system (1,2) with a broader power spectrum. Specifically, for identical coupling strength 
C = Co = 0.03, the index p = po (□) for system (1,2), while p = 0.3 > po (°) for system (3,4) although the frequency mismatch for system 
(3,4) is much larger. The effect of a Fourier band-pass filter applied to the system (1,2) while keeping C = 0.03 fixed is equivalent to an 
increase of the coupling strength of the system leading to a larger index pi > po (A) as also shown in Fig. 2(e). 



ity density function P(A?/'ii(i)) [Fig. 2(d)]. After band-pass 
filtering, A?/'i,i becomes smoother with less fluctuations, and 
the distribution P(A?/'i.i) exhibits a well pronounced peak. 
To quantify how the degree of synchronization changes with 
the width A/ of the band-pass filter, we calculate the synchro- 
nization index p [Fig. 2(e)]. We find that for very large values 
of the band width A/, the index p is the same as the value po 
obtained for the system (1,2) without any filtering, and that p 
remains unchanged for intermediate values of A/. However, 
for decreasing A/, the index p increases rapidly from the ex- 
pected value Po [Fig- 1(f). and 2(c)]. Such deviation to higher 
values of p > p^, while the coupling constant C in Eq. (1) 
remains fixed, indicates a spurious effect of synchronization 
due to the band-pass filter Thus, applying a band-pass filter 
with a too narrow band width when pre-processing empirical 
data may lead to overestimation of the phase synchronization 
(as defined by index p) between two empirical systems where 
the coupling strength is not known a-priori. 

Many physical and biological systems are influenced by ex- 



ternal noise, which can mask their intrinsic properties. Re- 
cent studies have shown that noise can bias the estimation of 
driver-response relationship in coupled nonlinear oscillators 
leading to change in synchronization measures [32]. Specifi- 
cally, external noise may weaken the detection of the coupling 
and reduce the synchronization between two coupled dynam- 
ical systems. To address this problem, we next test the effect 
of external noise on the degree of phase synchronization of 
the two coupled Rossler oscillators defined in Eq.(l). Adding 
uncorrected and unfiltered Gaussian noise rj to the output ob- 
servables xi and X2, while keeping the coupling constant C 
in Eq. (1) fixed, we find that the synchronization index p de- 
creases with increasing noise strength (j,^ (i.e., higher stan- 
dard deviation cr,, compared to the standard deviation a of 
the output signals xi and X2) [Fig. 3(a)]. The dependence of 
p on the value of the coupling constant G for different noise 
strength is shown in Fig. 3(b). We find that the transition to the 
state of maximum degree of synchronization (indicated by a 
horizontal plateau for p in Fig. 3(b)) occurs at decreasing val- 
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FIG. 2: Effects of band-pass filtering on synclironization. Time se- 
quence of the variables xi and X2 of system (1,2): (a) before and (b) 
after applying a band-pass Fourier filter with band width A/ — 0.01. 
After band-pass filtering the sequences xi and X2 are better aligned 
in time (with almost matching peaks), (c) Instantaneous phase differ- 
ence At/ji,!, and (d) the distribution P(AV'i,i) before (dashed line) 
and after (solid line) the Fourier band-pass filtering. After filtering, 
A?/)i_i is characterized by less fluctuations and a much narrower dis- 
tribution P{Atpi,i), indicating a stronger synchronization, although 
the coupling strength C — 0.03 remains constant, (e) Dependence 
of the index p on the band width 2nAf for fixed C — 0.03. A 
filter with a relatively broader band width (27rA/ > 1) leaves the 
synchronization index p practically unchanged, p = po, where po 
characterizes the synchronization between xi and X2 before filtering. 
Narrowing A/ leads to a sharp increase in p, which is an artifact of 
the Fourier filtering as the coupling C and all other parameters re- 
main unchanged, e.g, for A/ = 0.005, p = pi ~ 4po. 
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FIG. 3: Effect of external additive white noise on phase synchro- 
nization for system (1,2). (a) Dependence of the synchronization 
index p on the noise strength a,^ for fixed value of the coupling con- 
stant C. (b) Dependence of the synchronization index p on the cou- 
pling strength C for different levels of white noise which are defined 
through the standard deviation a^^ . 



ues of the coupling constant C for increasing noise strength 
ari- For very strong noise (ar^ — a = 8.3), the two Rossler 
oscillators in Eq.(l) appear not to be synchronized, character- 
ized by low values for the index p, even for very large values 
of the coupling constant C [Fig. 3(b)]. We note, that with in- 
creasing noise strength ct,, the position of the crossover to the 
plateau of maximum synchronization shifts to smaller values 
of C in Fig. 3(b), indicating that with increasing a^ the level of 
the plateau drops faster compared to the decline in the growth 
of p with increasing coupling C. 

To reduce the effect of noise in data analysis, a common 
approach is to apply a band-pass filter. In the case of the cou- 
pled Rossler oscillators defined in Eq.(l), we ask to what ex- 
tent a band-pass filter can reduce the effect of external noise 
while preserving the expected "true" phase synchronization 
as presented by po in Fig. 1(e). To answer this question, we 
first need to determine what are the limits to which spurious 
phase synchronization can be obtained purely as a result of 
band-pass filtering of two uncorrected and not coupled Gaus- 
sian noise signals. Our results for the synchronization index 
p obtained from multiple realizations of pairs of uncoupled 
white noise signals show that the synchronization index p 
can reach different maximum values Pmax^ indicated by ar- 
rows in Fig. 4(a), for different band width A/ — with de- 
creasing the band width pmax increases. The values of pmax 
provide an estimate of the maximum possible effect additive 
noise may have on the spurious "detection " of phase synchro- 
nization in coupled oscillators. Thus, empirical observations 
of synchronization index p > pmax may indicate presence 
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FIG. 4: Combined effects of external noise and Fourier band-pass 
filtering on the synchronization, (a) Cumulative distribution function 
F{p) = 1 - /J' P{p')dp' for the index p obtained from 100 differ- 
ent realizations of pairs of white noise signals without coupling. The 
length of the noise signals is inf [10^/27r]. Tails of the distributions 
for each band width indicate the maximum values of p one can obtain 
simply as a result of band-pass filtering when there is no synchroniza- 
tion between two white noise signals, (b) Synchronization index p 
obtained for system (1,2) with additive white noise as a function of 
the band width A/ for C — 0.03. While the effect of noise is grad- 
ually reduced by the Fourier band-pass filter with decreasing band 
width A/, there is an artificially increased synchronization (sharp 
increase in p) when 2nAf < 1, as also shown in Fig. 2(e). 



of a genuine phase synchronization between the outputs of 
two coupled oscillators, which is not an artifact of external 
noise. Our simulations show that the value of pmax does not 
change significantly with the length of the uncorrelated noise 
signals. In Fig. 4(b) we show how the synchronization in- 
dex p for system(l,2) depends on the strength of the added 
noise and on the width A/ of the band-pass filter. For very 
broad band width A/ the noise is not sufficiently filtered, and 
the synchronization between the two oscillators decreases (p 
decreases) with increasing noise strength cr^. With decreas- 
ing band width A/, i.e., applying a stronger filter, the effect 
of the noise is reduced, and correspondingly the index p in- 
creases — approaching the value po expected for the system 
(1,2) without noise. On the other hand, applying a filter with 
too narrow band width A/ leads to a spurious synchronization 
effects with p > po [Fig. 4(b)], following closely the depen- 



dence of p on A/ shown in Fig. 2(e) for a Rossler system 
without noise. 

In summary, our results indicate that phase synchronization 
between coupled nonlinear oscillators may strongly depend 
on the width of the power spectrum of these oscillators. Fur- 
ther, we find that while external noise can affect the degree of 
phase synchronization, band-pass filtering can reduce noise 
effects but can also lead to a spurious overestimation of the 
actual degree of phase synchronization in the system. This is 
of importance when analyzing empirical data in specific nar- 
row frequency ranges, for which the coupling strength may 
not be known a-priori. 
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